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Abstract 

A one-dimensional electromagnetic problem of Sommerfeld pre- 
cursor evolution, resulting from a finite rise-time signal excitation in 
a dispersive Lorentz medium is considered. The effect of the initial 
signal rate of growth as well as of the medium dumping on the pre- 
cursor shape and its magnitude is discussed. The analysis applied 
is based on an approach employing uniform asymptotic expansions. 
In addition, new approximate formulas are given for the location of 
the distant saddle points which affect local frequency and dumping of 
the precursor. The results obtained are illustrated numerically and 
compared with the results known from the literature. 

1 Introduction 

Fundamental investigations on EM signal propagation in the Lorentz model 
of a dispersive medium are due to Sommerfeld [Tj and Brillouin [21 E] • These 
authors revealed that in addition to the main signal propagating in the 
medium, two precursors are formed which precede the signal. The front of the 
fastest (Sommerfeld) precursor propagates in the medium with the velocity 
of light. The instantaneous oscillation frequencies of the precursors and their 
local dumping are directly related to the locations in the complex frequency 
plane of the corresponding saddle points in the integral representation of the 
signal. Those locations vary with space and time, and are governed by the 
saddle point equation (SPE) in Eq. 0, requiring that the phase function in 
the integrand in that representation be stationary. Analysis of the equation 
shows that there are two pairs of dominant saddle points, the distant and 
the near ones, responsible for the first (referred to as Sommerfeld) and the 
second (Brillouin) precursors, respectively. The saddle points in each pair 
are located symmetrically with respect to the frequency imaginary axis, the 
fact being related to the causality principle. 
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Since SPE seems to be not solvable in a closed form, attempts have been 
made to solve it in an approximate manner. Brillouin's approach [H] con- 
sisted in replacing the complex index of refraction in SPE by its expansion in 
powers of frequency, and then solving the simplified equation for frequency. 
As a result, a simple, approximate formula was found relating complex fre- 
quency to the distance and time coordinates. In the case of the Sommerfeld 
precursor, the applicability of this map was confined to the vicinity of the 
front of the precursor. Unfortunately, Brillouin, having at his disposal only 
a non-uniform asymptotic method, could not effectively describe the front 
evolution, as it corresponds to coalescence of two distant saddle points at 
infinity, the case not treatable with the method he used. This deficiency was 
removed by Bleistein and Handelsman j3] who developed a uniform asymp- 
totic approach that extended the validity of asymptotic considerations to the 
precursor front. 

Recently, the problem of asymptotic analysis of signal propagation in 
dispersive media has been reexamined by many authors. Extensive research 
in this field is due to Oughstun and Sherman j3] and Kelbert and Sazonov 
jH] . In a recent work, [7] Dvorak, Ziolkowski and Felsen offered a new hybrid 
approach combining both the asymptotic and FFT methods, the former being 
responsible for extreme parts of the signal frequency spectrum. 

In this work we reconsider the asymptotic model of signal propagation 
in a Lorentz medium, and concentrate on the Sommerfeld precursor. We 
assume that the Sommerfeld precursor is excited in the Lorentz medium by 
a sine modulated signal, with its envelope described by a hyperbolic tangent 
function. Such a representation provides a convenient model for signals with 
finite rise time. Unlike the Oughstun-Sherman study ((0, Sees. 4.3.4, 7.2.7 
and 7.3.7), where the signal envelope was described by an everywhere non- 
zero smooth function, our initial signal has the form of an abruptly switched 
modulated sine signal, i.e. it vanishes identically for t < and its envelope 
is non-zero for t > 0. At t = the derivative of the envelope suffers a step 
discontinuity. We study the influence of both the medium and the initial 
signal characteristics, including medium damping S and signal speed factor 
(3, on the evolution of the Sommerfeld precursor in the medium. In particular, 
we analyze how the speed factor affects the shape and the magnitude of the 
precursor excited by both the slow and the fast growing incident signals. We 
also obtain a simple approximation for the precursor damping factor. Finally, 
we provide a new approximation to the location of the saddle points, which 
is more accurate that those known in the literature. 

The results obtained here may appear to be useful e.g. in designs employ- 
ing fast Sommerfeld precursors as signals triggering the electronic devices 
designed to process the main signal. 
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2 The propagation problem and its exact so- 
lution 



We consider the ID problem of EM signal propagation in a Lorentz medium 
characterized by the complex index of refraction 



1/2 



n(0) 



1 



02 - 0^ + 2i50^ 



where 6^ = ATiNe'^/m is the plasma frequency of the medium, N, e and m 
are the number of electrons per unit volume, electron charge and its mass, 
respectively, 5 is a damping constant and 0o is the resonant frequency. 

In the plane z = of the Cartesian coordinate system {x, y, z}, the field 
component Ex{0,t), henceforth denoted by EqIJ:), is assumed to be given. 
It is described by a function, which has the form of a finite rise-time, sine 
modulated signal 



Eo{t) 



t < 

tanh(/3t) sin(0ct) t > 0. 



(2) 



The large positive coefficient (3 determines how rapidly the signal turns on, 
and 0c is a fixed carrier frequency. It is also assumed that no EM sources are 
present at z —>■ oo. In jS] a different signal was considered, also employing 
a tangent hyperbolic function, which, unlike (j21), did not vanish for times 
t < 0. 

In general, the problem consists in finding the field in the half-space 
z > and for time t > 0. The solution to this mixed, initial-boundary value 
problem for the Maxwell equations takes the form |Hj 



\ ria+oo 

E{z, t) = TT . 9(0 ; 0c) exp 



. C 



(3) 



where the constant a is greater than the abscissa of absolute convergence for 
EQ{t). The amplitude and phase functions g{0 ;P,0c) and \l/(0, 6*), respec- 
tively, are given by 



9i0;P,0c 




l[0 - 



2(3 



and 



- 0c (3 



m{0,e) = 0[n{0)-d]. 



i{0 + 0c) 
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(4) 
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The beta function is defined through the psi function as 

(6) 



(For the definition and properties of the psi function see jH], Sec. 6.3.) The 
beta function is related to the envelope of EQ{t) via the Fourier transforma- 
tion 

'ta„h^*e«*=iH(-|)-i. (7) 
Finally, 

.^f (8) 

is a dimensionless parameter that characterizes a space-time point {z, t). 

In this work we confine our interest to one component of the general solu- 
tion - the Sommerfeld precursor. A suitable approach to extract this partial 
field and study its dynamics is to evaluate the integral asymptotically as 
z oc. At the precursor onset, i.e. as 6^ — 1+, the distant saddle points in 
the complex plane meet at infinity to form a saddle point of infinite order. 
The valid asymptotic procedure that handles this case is a special instance 
of application of the general asymptotic theory developed by Bleistein and 
Handelsman and designed to uniformly evaluate integrals with nearby crit- 
ical points jini. (It can also be used in case of more than two coalescing 
critical points. The case of three critical points: a pole, a branch point and 
a saddle point was studied in jH].) 

Dynamics of the saddle points, essential in asymptotic considerations, is 
governed by the saddle point equation (SPE) 

n(0) + 0n'(0) - 6 = 0. (9) 

This equation results from the requirement that the phase © should be 
stationary. It has the form 6 = /(0). What we need is the inverse function 
= /~^(6'). Approximate solutions to the latter equation were obtained by 
Brillouin (Bj), Kelbert and Sazonov ([H]) and Oughstun and Sherman (j3]). 
In Sec. m we present another approximate solution, which is more accurate 
than the solutions known in the literature. 



3 Uniform asymptotic representation for the 
Sommerfeld precursor 

The phase function \E'(0,6') has the saddle point of infinite order at infinity 
(see pni, Ch. 9). In this case the classical asymptotic methods break down. 
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which imphes that they cannot be used to describe the precursor behavior at 
its front (i.e. for 9 —>■ 1^). Fortunately, the uniform approach, as proposed in 
|inj . can be effectively used. The term "uniform" means that the resulting 
asymptotic expansion is valid for any locations of the far saddle points in 
the complex plane, including the case where the points (symmetrical with 
respect to the imaginary axis) coalesce at infinity to create one saddle point 
of infinite order. As a consequence, the expansion holds for any 6' > 1, in 
particular at the precursor front. By using the methods described in [TO] . 
Sec. 9.5, we arrive at the following uniform asymptotic representation^ for 
the Sommerfeld precursor as 2; — > oo: 

E^(z,t) ~ -exp{-/ lm[^(0+,^)]} (10) 

X {Re[G'(0+;/5,0e)]Ji[-/ Re(^(0+,^))] 
+ Im[G(0+; /?, 0e)] J2H Re(^(0+, 6))]}, 



where 



G'(0+;/5,0e) 



Re[vI/(0+,^)] . , . . 

— — g[0 ,p,0cj, (11, 



\ ^00(0+, 

Ji(-) and J2(') are Bessel functions of the order 1 and 2, respectively, and 
I = z/c. 

In this paper we do not consider possible transition from the precursor to 
the main signal, which may happen for very high carrier frequency 0c- This 
occurs when the poles = ±0c — 2ikp, k = 0, 1, 2, . . . of the beta function 
are crossed during deformation of the original contour of integration in Q 
to the SDP through = 0^. 

A fundamental question that now arises is how the speed parameter /3 
in (j2I) and the damping parameter 6 affect the precursor dynamics in the 
medium. Below we try to answer this question. 

3.1 Dependence of the precursor on /3 

Let the Lorentz medium considered be described by Brillouin's parameters 



6 = 7201) X 10^S-\ 00 = 4.0 X 10^S-\ 6 = 0.28x10^^3-^ (12) 

and let us additionally choose 

0e = 2.0 X 10^^s-\ / = 5.0 X 10-^^s. (13) 

Assume further that f3 = 1.0 x lO^^s"^. The dynamic behavior of the Som- 
merfeld precursor field corresponding to this data, is shown in Fig. 1. 

"'^This representation is equivalent, but simplified in form, to that presented in J12'. In 
[H] the factor 1/4 in Eq.(4.5) should be replaced by 1/8. 
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Figure 1: Dynamic behavior of the Sommerfeld precursor in the Lorentz 
medium described by Brillouin's parameters. 0c = 2.0 x 10^^ s~^, I = 
5.0 X 10~^^s-i and (3 = 1.0 x lO^^s'^ 

If now (3 is increased, by one order or two, it is can be shown from (|TO|l that 
the oscillations in Fig. 1 are increased by the same factor, while the shape of 
the precursor is preserved. If, however, f5 takes much higher values, such as 
j3 = 1.0 X lO^^s"^ or more, the precursor shape distinctly changes (see Fig. 2), 
and the oscillation amplitudes virtually remain at the same, relatively high 
level as /3 further increases. 

This interesting behavior can be explained by studying properties of the 
function G(0"'"; P, 0c). The plot of G{0'^; P, 0c) as a function of P with fixed 
remaining arguments is presented in Fig. 3. 

The value 1.0001 of the variable 6 was chosen what corresponds to the 
close vicinity of the precursor front (cf. Fig. 1 and Fig. 2). It is seen from the 
plot that there are three characteristic regions of G{0~^] P,0c) variation. If 
for fixed 0~^ the parameter /3 is relatively small then the real part of function 
G{0'^] (3, 0c) is virtually zero and the essential contribution to the precursor is 
due to its imaginary part. This contribution increases in value with rising j3 
until about P = 4.0 x 10^ ''■s^^. This is the first region in which the precursor 
oscillation is described by the Bessel function J2. In the second, transitory 
region, the real part of G grows rapidly, and at about P = 2.0 x lO^^s^^ 
it settles down at a virtually constant level. The imaginary part reaches 
its maximum and then steadily decreases to zero. Here the contribution of 
the Bessel function Ji takes over. Finally, in the third region the real part 
of G remains nearly unchanged and the imaginary part vanishes. Now the 
contribution of Ji dominates and that of J2 is to be neglected. One can verify 
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Figure 2: Dynamic behavior of the Sommerfeld precursor in the Lorentz 
medium described by Brillouin's parameters. 0c = 2.0 x lO^^s"-*^; I — 
5.0 X lO-i^s"^ and /3 = 1.0 x lO^^s'^ 




/3 



l.xlO^" l.xlO^' l.xlO^" l.xlO^^ l.xlO 



Figure 3: Dependence of the real (solid line) and imaginary (dashed line) 
parts of the function G'(0+;/3,0c) on the parameter (3. Calculated for 6 — 
1.0001 and 0c = 2.0 x lO^^s'^ 
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that with increasing 6', the third region broadens, thus pushing the second 
region in the direction of smaller jS. 

Below we analyze the precursor behavior in the first and the third regions 
in more detail. To make the results as simple as possible, we shall employ 
Brillouin's approximation 

0+(e) ^ , ^ - 2i5. (14) 



The use of this approximation is justified in the vicinity of the precursor 
front. 

The case of relatively small P 

Here, by choosing sufficiently small 6 both arguments in functions ;B in Q 
can be made arbitrarily large. In this case the function G in (fTT|) may be 
simplified by substituting B for its asymptotic expression 

i3(w) = ^ + ^ + 0(M-^), M^±200, (15) 

valid also in some sectors centered around the rays argu = ±in/2. Then, 

G'(0+;/?,0e) 



\ ^,,{0+^9) '(0+^-02)2- ^ ^ 

If we expand this in fractional powers of 6' — 1 and retain the leading terms, 
we arrive at 

48v/2/3^0,(g - If' m0c{e - 1) 

G(0 ,/^,0e) p P . (17) 

Since 5\/6 — 1/6 is a small quantity, we can safely retain only the imaginary 
term. We see that G(0+;/3,0c) is proportional to {3. This fact confirms the 
observed behavior of G in the first region shown in Fig. 3. 

Let us now consider the envelope of the Sommerfeld precursor. It is 
obtained by replacing the Bessel functions in ()10|) by their envelopes. A 



good approximation for these envelopes is y2/(-iTx), where x denotes the 
argument of a Bessel function. Thus, by (fTO|l . the envelope of the precursor 
can be approximated by 

E^iz,t) ^ -exp{-/ lm[^(0+,^)]} 

^ ^ vr/ Re[^(0^^)] ^'^'^^^'^'^'''^^^ ^''^ 
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Figure 4: Sommerfeld precursor, its envelope and the envelope approximation 
at the precursor front. Calculated for /3 = 1.0 x lO^^s"^, / = 5.0 x lO^^^s and 
0c = 2.0 X lO^^s-i. 



Naturally, the notion of the envelope applies here to 6 greater than the first 
extreme of J2{x), occurring at x ~ 3.054. In our case it corresponds to 
e ^ 1.00009. 

If we are interested in the precursor front only, we can use (fT^ in (fTHjl . 
replace G{0'^; P, 0c) by RHS of (|T7jl . and expand the result in fractional pow- 
ers of 6' — 1. In this manner we arrive at the following approximation of the 
precursor dynamic behavior at its front: 



4 2V4/50c(e- 1)3/4(5 _ 65^2(0 -1))[1- 25/(^^-1)] 
^ ^) " ^ij^i • (19) 

In Fig. 4 an example of the precursor dynamics for (3 = 1.0 x 10^^s~^, its 
envelope as given by (fTHj) and the envelope approximation as given by (fT^ 
are shown. 

The slope of the envelope is given by 



dE^{z,t) _ 2"^(30,{h [3 - lAl5{e - 1)] + 6(5.^2(g - l)[18/^(^ - 1) - 5]} 



(20) 

For sufficiently small 6, the terms proportional to 6 can be neglected to yield 

rfE^(^,t)^, 3 21/4/50^ ^^^^ 



d9 b^'/^V^l {9 - 1)1/4 
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It is seen that the slope of the precursor envelope steadily decreases with 
growing 6. 

The case of large (5 

For finite 0+ and sufficiently large /3, the arguments in the functions i3 in Q 
can be made arbitrarily small. Then the asymptotic expansion is 

= - + ln2 + — + 0(^2), «^0. (22) 
u 12 

With its use in ffTTj) . the approximation for (^(0^; /3,0c), appropriate for the 
third region in Fig. 3, follows: 



G'(0+;/3,0e) 



Re(vl>(0+,g)) -0c [24/3^ + 7r^(0+^-0g)] 

\\ ^,,(0+,^) 12/32(0+^-02) • ^ ^ 



Proceeding as in the previous case, we find the following approximation for 
the envelope of the Sommerfeld precursor 

~s ^ [1 - 25l{e - 1)] [6^2(7 + 6)+ 4(g - 1) (96/3^ + 7,^01)] 0, 

192 21/4 53/2 (32^Q_ 1)3/4 

where minimal 9 (here, 9 ^ 1.00003) corresponds to the first extreme of 
Ji(x), occurring at x ~ 1.841. 

Fig. 5 shows the precursor dynamic behavior for f3 = 1.0 x lO^^s"^, its 
envelope and the approximation to the envelope as given by ()24j) . Notice that 
for any finite f3 there exists a 9, below which the assumption of large j3 and 
moderate 0^ is no longer valid. Therefore one should expect that with rising 
9, the precursor dynamics may pass through the stages described by the first 
or/and second regions in Fig. 3, before it reaches the stage characteristic of 
the third region. 



The slope of the precursor front is 

dE-^z, t) 9ti^0^{9 - 25 + 2(5/(3 + 29- h9^)] 



d9 768 21/453/2^ p2^Q _ 1)7/4 

? - 1)[1 - im{9 - l)](96/32 + 7r202) 



(25) 



+ 



192 21/453/2^7^/ (3^{9 - 1)7/4 



If terms proportional to the parameter 5 are neglected, this expression re- 
duces to 

dE^{z, t) 0,[b^7c\9 - 25) + A{9 - 1){96(3' + 7c'0l)] 

— . (26) 



d9 768 21/453/2^^ p2(^g _ 1)7/4 
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Figure 5: Sommerfeld precursor, its envelope and the envelope approximation 
at the precursor front. Calculated for j3 = 1.0 x 10^^ s~^, I = 5.0 x lO^^^s and 
0c = 2.0 X 10^^s-\ 



As before, the envelope slope decreases with 6. For the parameters used here 
the rate of the precursor growth is about 16 times higher than in the previous 
case. 

If we formally let f3 tend to infinity, we obtain from ()25|1 

dE'iz, t) ^ 23/40,[l - 1075/ {e - 1)] ^27) 



de 63/2v%/ {e - 1)3/4 

which corresponds to the unit-step function envelope in the initial signal. 

Finally, let us consider the precursor behavior at the first oscillation, 
provided (3 is finite. With the use of (dH), (fTTj) and power expansions of the 
Bessel and exponential functions, we find for ~ 1 (see Fig. 6) 

E%z,t)^(3l0,{e-lf[l-'^-^Y (28) 

This implies that 

, r.?/ X , dE^(z,t) , d^E^(z,t) , 
lim E^(z,t) = lim = 0, lim ^ 0. 

On the other hand, we obtain 

, ^ , , , dE^it) , d^E^it) , 

limEo t = lim — = 0, lim — ^ 0. 

t^o '^^ ^ t^o dt t^o dt^ 

It then follows that the Sommerfeld precursor has the same smoothness prop- 
erties about 6' = 1 as the initial signal Q has about t = 0. This confirms 
the more general results based on Green's function approach ^| and [TIj . 
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Figure 6: Onset of the Sommerfeld precursor calculated for = 1.0 x lO^^s^"*^, 
/ = 5.0 X lO^-'^^s and 0c = 2.0 x lO^'^^s""'^. The dashed line approximates the 
precursor growth near 6=1^ (Eq. W^) . 



3.1.1 Special case 

Assume that /3 oo and 6^ ~ 1+. Then 51(0+; /5,0c) -20^^2(6* - 1)6" 
and (frnjl reduces to 



E'{z,t) ^ -^—^ -J,{lbp{6-1)). 

It is readily seen that this result fully agrees with the representation 

i?^(^,t)^^^Ji(2v/^^) 

obtained by Sommerfeld on the grounds of integral considerations and valid 
for the initial signal described by the Heaviside unit step function (,3j, Eq. (33)). 
Here, we have employed Sommerfeld's notation: 

z ^ b'^z 2tt 

t = t- - ^ = r = — . 

c 2c 0c 

3.1.2 A comment on the form of asymptotic representation of the 
precursor 

The form of the asymptotic representation of the Sommerfeld precursor de- 
pends on the way Eo{t) behaves at t ~ O"*". First we note that, (^15J, Eqs. 
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(7.128)), the asymptotic behavior of EQ{t) just after it is turned on: 




(29) 



imphes the following asymptotic behavior of g at infinity: 



9(0;/3,fc)~|;(^) 



r+1 



as 







oo. 



(30) 



If the RHS of (jHn|l is used in the basic integral formula describing the signal 
evolution in a dispersive medium, it appears ([IS], Eqs. (7.144)) that in the 
vicinity of the front, i.e. for 6 1+, the Sommerfeld precursor dynamics is 
described by Jr[b^j2{9 — l)z/c]. If (3 oo, one has r = 1 (which follows from 
expanding the sine function alone), and the precursor dynamics is described 
in terms of Ji. If P is finite, r = 2 (see below (|1H1)), and now J2 describes 
the precursor development. The orders of the Bessel functions appearing in 
the uniform asymptotic representation of the precursor are determined from 
a similar criterion, relating these orders to the behavior of the integrand at 
infinity. In our case the smallest value of r is 1, which implies the presence 
of the functions Ji and J2 in the asymptotic formula (fTUj) . 

With increasing 6, the real value of 0"'" in the arguments of beta func- 
tions in (jlj decreases. This results in moving the boundaries between the 
three characteristic regions of G{0'^; f3, 0c) variation with f3 in the direction 
of smaller values of p. As a consequence, Ji describes the precursor behav- 
ior not only at /3 — 00, but also at finite, sufficiently large values of this 
parameter (see Fig. 3). 

3.2 Decay of the precursor 

By expanding the phase function \E' in terms of powers of 6 we obtain 



Since the imaginary part of the complex frequency 0+ is small compared to 
the real part of this frequency, the imaginary part of the phase function 
can be approximated with the second term in this expansion. 

For 9 small, |0+| >> 6, and we neglect 6 in the saddle point equation 
to arrive at 



^(0+,^) = 




WiS^). (31) 



6^0+^ - 02 0+- 




(32) 



(0+^ - 02)2 + 0+^ - 02 0+^ 
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Figure 7: Comparison of Im "^{0^,6) (solid line) with the approximations 
given in (Eq Ul]) (dash-dash line) and fEqU3) (dash-dot line), respectively. 

For not too large, the square root at the RHS in this equation, and conse- 
quently the second term at its LHS, as well as the square roots in (j3T| . can 
be approximated by 1. It follows then that the first term in (jS21)j appearing 
also as a factor in the second term of (jSH), equals approximately 6 — 1. Thus 
by (EH), 

lm^(0+,e) ^ 1). (33) 

The function Im[\E'(0+, 9)], and its approximations given in ()31|) and ()33|) are 
shown in Fig. 7. 

It follows that major cause of the precursor decay is its exponential damp- 
ing. For low and moderate 9, the exponent determining the rate of decay is 
approximately equal to —16{9 — 1). The contribution to precursor decrease 
due to natural lowering of the oscillation amplitudes in the Bessel functions 
is of smaller importance. 

4 The location of the saddle point 0^(0) 

As seen from fllOj) . evolution of the precursor depends on the function 0'^{9) 
which is a solution to SPE, describing the location of the distant saddle point 
of (0) in the right complex half-plane. Since this function does not seem 
to be expressible in a closed form, approximate formulas for it were found 
in the literature (Brillouin 3J, Oughstun and Sherman [S], and Kelbert and 
Sazonov jO]). In this section we offer another approximate solution to (jH)), 
which is more accurate than the previous ones. 
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We begin with approximate solution to SPE found in [T^ which provides 
good accuracy for large and medium 6, but fails for 6 close to 1. First let us 
notice that 0^{9) can be written down in the form 



where 



i-n\e) 



62 



The approximation obtained in [I6j results from substituting for n{6) 



n g 



1 (e 



2c] -g^ 



\ 4a \g 

which leads to approximation of C(^)) to be denoted by (2- In fl36|) . 
1 



(34) 
(35) 

(36) 



9 



2-1/3 



{u + Vu^ - t;3)2/3 + 



u = 2c^ - 72ace + 27a^^ 



-2c, 



(37) 
(38) 



and 



V = 2^/^(0^ - 12ae). (39) 
The coefficients a, c and e are constant for a given medium and are equal to 



a 



i5J0i - s 



(40) 



i6 



2J0I - 52 62 



8(02 _ 52)3/2 
0l-5^-r5^0l-5^y (41) 



and 



0? - 52 _ ^5^/02 _ 52 



(42) 



We shall denote this approximation by 0'^{6X2)- 

For 9 1+, i.e. in the vicinity of the precursor front, we find here another 
approximation. By expressing SPE in terms of the variable C, we obtain 



On -I 



(43) 
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where w 



— 6"^. If ^ — > 1+ then ( and this equation can be 



approximated by 



e-1 9 



(44) 



The solution to ()44j) . relevant to our problem can be found by means of the 
Mathematica computer program: 



2w 



u 



+ 



2m2 - 3s 



4m 



(45) 



where 



V 

w 

I = 6'^ + ew, u 



J = b'e' -8{e-i)w, h = -w{e-i), 



(46) 



J 



p 



p 



V 4 662' 6^2 \^22/3p 24/3 J ' 

(r + Vr2 -4^3)^/3^ z = b^[f - 4 12b^h + 4^ 12w'(^ - 1)^] 



r = b''[2j' - - 12'{e - lyiSjw' - 1262/2 _ i262^^^^^)]. 

By substituting ( for in (jSl)), a new approximation of the distant saddle 
points location is obtained which is valid for 6 close to unity. We denote it 

by 0+(^,Ci). 

The two approximations, 0~^{9,Ci) and 0~^{9,(2), can now be combined 
into one formula that provides smooth transition from one approximation to 
the other. For Brillouin's choice of the medium parameters we choose the 
transition value of 6 to be 1.3. Then the joint approximation can be written 
down as 



05d(^) = 0+(^;Ci) 



H{i.3-e 



sign(e - 1.3) k{9 - 1.3) 



Hi 



2 

1.3) 



<0) 

sign(^ - 1.3) /€(^ - 1.3) 



(47) 



n{0) 



where H{s) is a unit-step function, 

k{s) =77(5 + 0.05)77(5 -0.05), 

and 



exp(-l/s2), s > 0, 



0, 



s < 0. 



(48) 
(49) 
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Since k{9 — 1.3) is zero outside the interval 1.25 < 6 < 1.35, the approxi- 
mation 0sd{O) equals 0+(^,Ci) if ^ < 1-25, and 0+(^,C2) if ^ > 1-35. The 
factors in square brackets in (j47jl provide smooth transition between the two 
approximations inside the interval. At ^ = 1.3, they are understood in a 
sense of limits (from the left or from the right), and thus are equal to 1/2. A 
different choice of media parameters may require the numerical parameters 
in dHI) and (gHl) to be modified. 



Re (w) 




Figure 8: Real and imaginary parts of various approximate solutions to 
(Eq\^. Legend: E - exact result (solid line), B - Brillouin's approxima- 
tion, OS - Ougstun- Sherman's approximation, KS - Kelbert-Sazonov's ap- 
proximation, N - the result based on (Eq\47\), and ^2 ^ the approximations 
described by0^{6,^i) and 0^^ {9 , (,2) , respectively. 

The approximation ()47|1 is shown in Fig.8, together with the approxima- 
tions obtained by Brillouin, Kelbert and Sazonov, Oughstun and Sherman, 
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and with the partial approximations 0^(6'; Ci) and 0^(6*; C2)- This approxi- 
mation is not as simple as the first three ones, but may prove to be useful 
when a higher accuracy is required. Its maximal deviation from the solution 
of found numerically slightly exceeds 2 percent. 

Note that the approximation 0'^{6; (2) indicates an interesting symmetry 
between locations of the distant and near saddle points (see Similar 
symmetry also follows from the Kelbert-Sazonov approximation. 
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